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ABSTRACT 

We present NuSTAR observations of the bright Seyfert 2 galaxy NGC 2110 obtained 
in 2012, when the source was at the highest flux level ever observed, and in 2013, 
when the source was at a more typical flux level. We include archival observations 
from other X-ray satellites, namely XMM-Newton , Suzaku , BeppoSAX , Chandra and 
Swift. Simultaneous NuSTAR and Swift broad band spectra (in the 3-80 keV range) 
indicate a cutoff energy E c > 210 keV, with no detectable contribution from Compton 
reflection. NGC 2110 is one of the very few sources where no evidence for distant 
Compton thick scattering is found and, by using temporal information collected over 
more than a decade, we investigate variations of the iron Kq line on time scales of years. 
The Fe Ka line is likely the sum of two components: one constant (originating from 
distant Compton-thick material) and the other one variable and linearly correlated 
with the source flux (possibly arising from Compton-thin material much closer to the 
black hole). 

Key words: Galaxies: active - Galaxies: Seyfert - Galaxies: accretion - Individual: 
NGC 2110 


1 INTRODUCTION 

The X-ray spectra of Seyfert 2 galaxies offer a unique op- 
portunity to probe the circumnuclear environment of Ac- 
tive Galactic Nuclei (AGN). The obscuration of the pri- 
mary radiation by matter with column densities typically 
in the 10 22 — 10 24 cm -2 range allow a study of the cold 
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and ionised reflectors that cannot be observed, due to dilu- 
tion effects, in unobscured sources. Typical X-ray features 
of the cold circumnuclear material include an intense Fe 
Kq line at 6.4 keV due to fluorescence e mission a nd a re- 
flection continuum pea king at ~ 30 keV llMatt et alJll99ll : 
iGeorge fc Fabiarilll99ll ). 

The primary continuum is thought to arise from the in- 
nermost regions surrounding the central black hole, in a hot 
corona above the accretion disc. It is observed as a power-law 
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spectrum extending to energies determined by t he electro n 
temperature in the hot corona (iRvbicki fc Lightmanll 19791) . 
The power-law index is a function of the plasma temperature 
T and optical depth r. 

NGC 2110, at a redshift z = 0.008, is one of the bright- 
est Seyfert galaxies in the hard X-ray (>10 keV) band and 
it shows a promine nt Fe Kq line ac c ompanied by variabl e 
intrinsic emission (iMushotzkvl 1 19821 ; lHavashi et al.l 1 1996h . 
It has been intensively studied by most X-ray observato- 
ries, and it has shown very inte resting and peculiar char- 
acteristics. iMalaguti et al.l (Il999l l reported a photon index 
of T = 1.86 analysing BeppoSAX PDS data above 13 keV 
(consistent with typical values found in Seyfert 1 sources: 
iNandra fe Poundall994fl . However, when 2-10 keV data are 
considered, the photon index becomes flatter ( T = 1.67) 
and the authors suggested the presence of obscuring material 
with a column density Ah = 4.1+q'| x 10 22 cm -2 along the 
line of sight and an additional absorber with Ah ~3x 10 23 
cnh 2 , partially covering the nuclear source. This scenario 
was conf irmed by XMM-Newton , Chandra and Suzaku obser- 
vations fevans et al.l 120071 : iReeves et al.l 120061 : iRivers et al.l 
1201411 . 

NGC 2110 is one of the very few Seyfert galaxies that, 
despite the intense iron K a emission line at 6.4 keV, does not 
show any Compton reflection from circumnuclear material: 
values of R ^ 0.17 and R T 0.1 were found with BeppoSAX 
and S uzaku , respectively (IMalaguti et al.lll999l : IRivers et al.l 
l2014l l. If the line emitting material is Compton-thick (Ah > 
10 24 cm -2 ), the iron Ka emission would be accompanied by 
a Compton reflection compon ent above 10 keV dMatt et alJ 
Il99ll : ICeorge fc Fabianlll99ll l. The inferred upper limits on 
the reflection fraction R in this source suggest that the iron 
Kq line is not produced by distant, Compton-thick material 
but is instead emitted by Compton-thin matter, such as in 
the case of NGC 7213 llBianchi et al.ll2003l . 120081 ). 

Recently, NuSTAR observed NGC 2110 in an extremely 
bright state. The lack of reflection components in NGC 
2110 arising from the accretion disc or from the putative 
torus above 10 keV makes this source a perfect candidate 
for measuring a high-energy cutoff. We present a detailed 
multi-epoch X-ray study of NGC 2110 including two recent 
NuSTAR observations, obtained in 2012 and 2013, with the 
dual aims of studying the behaviour of the Fe Kq line with 
respect to the highly variable intrinsic continuum and of 
constraining the high energy cutoff in this source. 

2 OBSERVATIONS AND DATA REDUCTION 

2.1 NuSTAR 

NuSTAR (Harrison et al. 2013) observed NGC 2110 with 
its two coaligned X-ray telescopes with corresponding fo- 
cal planes: Focal Plane Module A (FPMA) and B (FPMB) 
starting on 2012 October 5 and 2013 February 14 for a to- 
tal of ~ 32 ks and ~ 26 ks of elapsed time, respectively. 
The Level 1 data products were processed with the NuS- 
TAR Data Analysis Software (NuSTARDAS) package (v. 
1.3.0). Cleaned event hies (level 2 data products) were pro- 
duced and calibrated using standard filtering criteria with 
the NUPIPELINE task and the latest calibration Hies avail- 
able in the NuSTAR calibration database (CALDB). Both 


extraction radii for the source and background spectra were 
1.5 arcmin. After this process, the net exposure times for the 
two observations were 15 ks and 12 ks, with 3-80 keV count 
rates of 6.53 ± 0.02 and 4.50 ± 0.02 cts s -1 for FPMA, and 
6.32 ±0.02 and 4.25 ±0.02 cts s -1 for FPMB. The two pairs 
of NuSTAR spectra were binned in order to over-sample the 
instrumental resolution by at least a factor of 2.5 and to 
have a Signal-to-Noise Ratio (SNR) greater than 5 a in each 
spectral channel. 


2.2 Suzaku 

NGC 2110 was observed by Suzaku starting on 2005 Septem- 
ber 16 (OBSID 100024010) and seven years later, starting on 
2012 August 31 (OBSID 707034010). Data were taken from 
the X-ray Imaging Spectrometer (XIS) and the Hard X-ray 
Detector (HXD). We reprocessed the event Hies with the 
latest calibration files available (2014-02-03) using FTOOLS 
6.14 and SUZAKU software Version 21, adopting standard 
procedures and recommended screening criteria. The source 
extraction radius was 1.8 arcmin. Background spectra were 
extracted from source-free regions of 1.8 arcmin radius. Re- 
sponse matrices and ancillary response files were generated 
using XISRMFGEN and xisarfgen. The 0.5-10 keV spectra 
extracted from the front-illuminated XIS0 and XIS3 were 
co-added using the ftool ADDASCASPEC, for net exposure 
times of 102 ks and 103 ks for the two data sets. Spectra 
were binned in order to over-sample the instrumental reso- 
lution by at least a factor of three and to have no less than 
30 counts in each background-subtracted spectral channel. 
This allows the applicability of \ 2 statistics. 

We reduced the HXD PIN data using the AEPIPELINE 
reprocessing tool, and for background determination we 
downloaded and utilized the tuned non- X-ray background 
(NXB) provided by the HXD team. We extracted source 
and background spectra using the same good time intervals. 


2.3 XMM-Newton 

NGC 2110 was observed by XMM-Newton (I Jansen et alJ 
l200lll for ~60 ks, starting on 2003 March 5 (OB- 
SID 01456701 01) w ith the EPIC CC D cameras, t he Pn 
dStruder et aLll200ll l and the two MOS dTurner et alJliooil l. 
operated in large window and thin filter mode. Data from 
the MOS detectors are not included in our analysis due to 
the lower statistics of the spectra. The extraction radii and 
the optimal time cuts for fl aring part icle background were 
computed with SAS 13 (iCabriel et al.ll2004T ) via an iterative 
process which leads to a maxi mization of the S NR, s imilar 
to the approach described in IPiconcelli et al.t ( 2004 1. The 
resulting optimal extraction radius was 38 arcsec and the 
background spectra were extracted from source-free circular 
regions with a radius of about 50 arcsec. 

After this process, the net exposure time was 44 ks for 
the EPIC-P 11 . Spectra were binned in order to over-sample 
the instrumental resolution by at least a factor of three 
and to have no less than 30 counts in each background- 
subtracted spectral channel. 


The hard X-ray emission of NGC 2110 observed by NuSTAR 3 


Obs. ID 

Date 

Texp (ks) 

35459001 

2006-03-25 

8.3 

35459002 

2006-04-08 

9.2 

35459003 

2006-04-15 

2.3 

35459004 

2008-08-31 

2.2 

35459005 

2009-10-12 

3.5 

80364001* 

2012-10-05 

7.1 

35459006 

2013-02-03 

14.1 

80364002* 

2013-02-15 

0.9 

35459007 

2013-03-09 

2.7 


Table 1. Observation log for the Swift monitoring of NGC 2110. 
Observations IDs, dates and net exposure times (after filtering 
and correction for photon pile-up) are reported. Asterisks indicate 
observations simultaneous with NuSTAR. 


2.4 Swift 

There were 9 Swift/XRT observations of NGC 2110, in 2006, 
2008, 2009, 2012 and 2013. The first 8 observations were car- 
ried out with the photon counting (PC) mode, while the last 
observation was in window timing (WT) mode (Table [p. 
ObsIDs 80364001 and 80364002 are simultaneous with NuS- 
TAR. We reprocessed all the datasets to generate cleaned 
event files with the xrtpipeline script included in the hea- 
SOFT version 6.13. The PC mode observations are all af- 
fected by pile-up, while the final observation, obtained in 
WT mode, does not have pile-up issues. To correct the pile- 
up, we extracted the source spectrum within a 72 arcsec 
radius circular region, excluding the central 10 arcsec radius 
aperture for the PC mode datasets. Obs. IDs, dates and net 
exposure times are listed in Table |T] 


2.5 Chandra 

NGC 2110 was observed by Chandra on 2000-04-22 with the 
ACIS-S camera (Obs. ID 883), on 2001-12-19 with three con- 
secutive HETGs pointings (Obs. ID 3143, 3417 and 3418) 
and on 2003-03-05 with the HETGs (Obs. I D 4377). A ll 
t hese observations are discussed in detail in I Evans et al.l 
(120071) where the four HETG spectra were merged. Since 
pile-up was found in the ACIS-S observation, we do not use 
this spectrum in our analysis. Data were reduced using CIAO 

4.5 and the latest CALDB files, we merged Obs. IDs 3143, 
3417 and 3418 spectra using the add_grating_ORDER 
and ADD —GRATING— SPECTRA tools. The resulting exposure 
times are 105 ks and 95 ks for the 2001 and 2003 HETGs 
data, respectively. Data were binned to have no less than 30 
counts in each spectral channel. 


2.6 BeppoSAX 

The source was observed by BeppoSAX on 1997-10-12 with 
the MECS for a net exposure time of 76 ks. Reduced data 
were downloaded from the BeppoSAX online data archive 0- 


1 available at: 

http : //www. asdc . asi . it/mmia/index . php?mission=saxnf i 


Confidence contours: Chi-Squared 



Figure 1. Contour plot between reflection fraction R and photon 
index 1 for the high flux NuSTAR observation in 2012, when 
only 10—80 keV data are considered. The solid black, red and 
green lines correspond to 68%, 90% and 99% confidence levels, 
respectively. 

3 SPECTRAL ANALYSIS 

We first study the simultaneous NuSTAR and Swift data to 
probe the primary radiation from the AGN and the proper- 
ties of the hot corona. We then perform a multi-epoch phe- 
nomenological X-ray analysis to study the behavior of the 
Fe Ka emission line in response to variability of the nuclear 
continuum. 


In previous work, broad band analyses of some data 
sets revealed the presence of extra-nucl ear emission in the 


softer (E<1 keV) 

Dart of the spectra (lEvans et al. 

200(| 

Reeves et al. 20061; Evans et al. 2007; Rivers et al. 2014 

), the 


analysis of this component is beyond the scope of this work. 
We focus our analysis on the 3-79 keV band where the con- 
tribution from soft diffuse emission is negligible. 


The adopted cosmological parameters are Ho = 70 km 
s -1 Mpc -1 , Ba = 0.73 and Q, m — 0.27, i.e. the default ones 
in XSPEC 12.8.1 (lArnaudlll996l ). Errors correspond to the 
90% confidence level for one interesting parameter (Ay 2 = 
2.7), if not stated otherwise. 

3.1 NuSTAR 3—80 keV spectral analysis 

The X-ray spectra of NGC 2110 have been widely studied 
in the past few years and several analyses have shown the 
presence of a T ~ 1.7 power law continuum partially covered 
by several layers of absorbing material with column densities 
in the range 2-6 xlO 22 cm -2 llEvans et al.ll2007l : lRivers et alJ 
l2014ll . Above 10 keV the effect of this absorbing material 
can be neglected and we therefore start our spectral analysis 
fitting the NuSTAR spectra in the 10-80 keV range, to have 
a direct measurement of the photon index of the primary 
power law. We will test a posteriori if this assumption is 
correct. 

Our model is composed of a power law, with the only 
free parameters in the fit being the photon index, the 
normalization of the power law and the cross-calibration 
factors between the FPMs. When the model is applied 
to the NuSTAR 2012 and 2013 spectra, the fit is good 
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NuSTAR FPMA/B and Swift XRT best Tit 



Figure 2. 3-80 keV simultaneous NuSTAR-\- Swift best fit. NuS- 
TAR FPMA and FPMB spectra are in black and red for the 2012 
observation, green and blue for the 2013 one. Swift XRT 2012 
data are in light blue and 2013 spectra in magenta. Residuals are 
shown in the bottom panel, when a model composed of an ab- 
sorbed cutoff power law and a Gaussian line at 6.4 keV is applied 
to our data set. No strong residuals are present accross the whole 
energy band. 


(x 2 /dof=568/601=0.95) and no strong residuals are present 
across the whole energy band. The best fit photon index 
isT = 1.691 ±0.015 and the cross-calibration factors are 
Kfll = 1.017 ± 0.015 and Kl°4| = 1.000 ± 0.020. If we 
leave the photon indices of the 2012 and 2013 observations 
free to vary independently no improvement in the fit is found 
and they are both consistent with the inferred single best fit 
value. A neutral absorber along the line of sight does not 
improve the fit (Ay 2 = 2 for one additional variable param- 
eter) and only an upper limit Ah < 8 x 10 22 cm -2 is found: 
this confirms that the circumnuclear absorbers in this object 
do not affect the analysis of the 10-80 keV NuSTAR spectra. 

In the past, high energy observations of NGC 2110 only 
revealed upper limits for the fraction R of the Compton re- 
flection of the primary continuum by cold , neutral material 


llMalaguti et al.l Il999l : iRivers et al. 2014li which is usually 
found in Seyfert galaxies ( Dadinal 20081 : Ricci et~ahl 1201 if) . 
The addition of such a component in our fit, modeled with 
pexrav (IMagdziarz fc Zdziarskilll995f) with T fixed to that 
of the primary continuum and the inclination angle to 60 de- 
grees, does not lead to an improvement of the fit (Ay 2 = 2 
for two additional variable parameters) and only upper lim- 
its -R2012 < 0.15 and R2013 < 0.14 in the reflection compo- 
nents are found for the 2012 and 2013 observations, respec- 
tively. The contour plot between R 2012 and T is shown in 
Fig. [T] Other parameters do not differ from best the fit val- 
ues presented above. When we add a cutoff power law in our 
fit (model CUTOFFPL in xspec) no significant improvement 
is found (Ax 2 = 2 for one additional variable parameter): 
we find a best fit value of T = 1.647 ^q' 353 and a lower limit 
E c > 230 keV for the cutoff energy. 




Figure 3. Top panel: Contour plot between photon index F and 
cutoff energy E c , when simultaneous 3-80 keV Swift.+NuSTAR 
data are considered. The solid black, red and green lines corre- 
spond to 68%, 90% and 99% confidence levels, respectively. High- 
energy cutoff values are limited to ^500 keV because it is the 
maximum value allowed in the CUTOFFPL model. Bottom panel: 
Contour plot between photon index F and absorbing column den- 
sity along the line of sight. 


3.2 NuSTAR+Swift 3—80 keV spectral analysis 

We then consider 2012 and 2013 NuSTAR data down to 
3 keV and introduce the simultaneous Swift-XRT spectra, 
with net exposure times of 7.1 ks and 0.9 ks, respectively. 

The model is composed of an absorbed power law 
with a high-energy cutoff and two Gaussian lines, to re- 
produce the Fe Ka and K/3 emission at 6.4 keV and 7.056 
keV, respectively. We then convolved the baseline model 
with a Galactic colu mn density Ah = 1.62 x 10 21 cm -2 
llKalberla et al.ll2005l), modeled with tb abs in xspec, us- 
ing s olar abundances (IWilms et al.ll 2000 l l and cross-sections 
from IVerner et ahl d 19961 ). In XSPEC the model reads as fol- 
lows: TBABS*ZWABS*(CUTOFFPL+ZGAUSS-|-ZGAUSS). XRT- 
FPMA cross-calibration factors are introduced as a variable 
parameters. The fit is good (x 2 /dof=876/881=0.99) and no 
strong residuals are seen (Fig. 0. 

The inclusion of Compton reflection in the fit leads to 
an insignificant improvement of the fit (Ay 2 = 2 for two 
additional variable parameters) and its contribution to the 
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Instrument 

Date 

A h 

r 

Fe Ka En. 

C 7 

EW 

Fkc, 

Fk/3 

F 3 - 10 keV 

BeppoSAX 

1997-10-12 

4.3 ±0.9 

1.74 ±0.09 

5 40+O.O6 
0.09 

< 280 

194^50 

0 0+3. 0 
"•^-2.3 

< 1.3 

2.77 ±0.05 

Chandra 

2001-12-19 

4.0 ± 1.8 

1 07+0.30 
1 ' u -0.25 

6.400 ± 0.008 

10+ 14 

lo -10 

90^25 

5 4+ 1 - 8 

< 2.8 

3.84 ± 0.07 

Chandra 

2003-03-05 

< 4.5 

i qc+0.48 

x.^o_o 33 

6.391 ±0.016 

30+ 31 

ou -16 

135±°° 

c c+2.5 
°-°-2.0 

< 2.0 

2.80 ±0.07 

XMM 

2003-03-05 

3.9 ±0.4 

1.57 ±0.05 

6.42 ± 0.01 

62 ± 14 

145 ± 15 

5.0 ±0.5 

0.8 ±0.3 

2.26 ±0.03 

Suzaku 

2005-09-16 

3.8 ±0.2 

1.63 ±0.02 

6.40 ±0.01 

50 ± 15 

55 ± 5 

8.4 ±0.8 

0.6 ±0.5 

9.90 ±0.03 

Suzaku 

2012-08-31 

4.5 ±0.2 

1.63 ±0.02 

6.39 ±0.01 

< 55 

50 ±7 

9.7 ±0.9 

< 1.0 

11.8 ±0.1 

NuSTAR 

2012-10-05 

4.0 ±0.4 

1.64 ±0.03 

6.33 ±0.07 

< 192 

35 ± 10 

9.5 ±3.0 

< 2.3 

17.1 ± 0.2 

NuSTAR 

2013-02-14 

4.0 ±0.7 

1.64 ±0.05 

6.45 ± 0.07 

1 7C+2 00 
- uo -140 

90 ±% 

I6.0±|;g 

< 3.4 

11.7 ± 0.2 





Swift best fit parameters 





X 2 /d.o.f. 

Date 

A h 

r 

Fe Ka En. 

a 

EW 

Fka 

Fkp 

7*3- 10 keV 

59/61 

2006-03-25 

3.5 ±2.5 

1.35 ±0.45 

6.50 ±0.15 

H5t 8 8°5 

190 ± 115 

25.0 ± 15.0 

< 13 

8.9 ±0.5 

83/85 

2006-04-08 

3.0 ±2.0 

1.25 ±0.35 

D.DU_o os 

< 250 

100 +8 ° 

luu -64 

14.7±i° 3 2 

< 7 

9.8 ±0.5 

15/24 

2006-04-15 

4.0 ±2.3 

1.7* 

6.7 ±0.5 

< 350 

450~*” 665 

^ ou _400 

5o.9tS; S 5 

< 34 

8.1 ±0.9 

17/15 

2008-08-31 

7.5 ±3.0 

1.7* 

6.4* 

60* 

< 95 

< 9.0 

< 44 

5.9 ±0.9 

37/35 

2009-10-12 

5.5 ±3.5 

1 40 + O.7O 
— 0.51 

6.2 ±0.1 

< 350 

170 -1-120 
1 ‘ u — 120 

35 ±25 

< 38 

11.4 ±0.8 

106/94 

2012-10-05 

5.9 ±2.5 

1 7+O.3 

6.4* 

60* 

< 35 

< 30.7 

< 15 

17.3 ±0.4 

170/159 

2013-02-03 

6.3 ± 1.9 

1 '7C+0.25 
'°-0.25 

6.4* 

60* 

< 50 

< 14 

< 13 

16.0 ±0.9 

9/8 

2013-02-15 

7.0 ±3.5 

1.7* 

6.4* 

60* 

< 250 

< 53 

< 30 

12.5 ± 1.6 

81/56 

2013-03-09 

4.3 ± 1.2 

1 . 7 * 

6.2 ±0.1 

60* 

160 ± 100 

22.7 ± 12.5 

< 10 

8.0 ±0.8 


Table 2. Best fit parameters when the 3-10 keV phenomenological fit is applied to the data. Column densities are in 10 22 cm -2 units, 
energy centroids are in keV units, equivalent widths (EW) and widths (rr) of the Fe Ka line in eV units, line fluxes are in 10” 5 ph cm -2 
s' 1 units and 3-10 keV observed fluxes are in Hr 11 erg cm -2 s — 1 units. 


total 3-80 keV flux is Ffff 2 < 1.5% and F 2 ff 3 < 2.5%. Best 
fit cross-calibration factors are ^xrt-fpma = 0-98 i 0.03 
and ^xrt— fpma = 0.95 ± 0.09. We measure a lower limit 
for the high energy cutoff E c > 210 keV and in Fig. [3] the 
contour plots of E c vs T as well as Ah vs T are shown. Best 
fit values do not significantly differ from the ones discussed 
in Sect. 13.11 

We next use a more physical model (compTT in XSPEC; 
iTitarchu^ ll994T ) to measure the coronal temperature kT e 
and optical depth r. In this model the soft seed photon spec- 
trum is a Wien law; we fixed the disc temperature to 30 eV, 
appropriate for Mbh ~ 10 s (the black hole mass in 
NGC 2110 is estimated to be Mbh ~ 2 x 10 8 Mq via the 
Mbh-u relation; iMoran et alj l2007l l . In the case of a slab 
geometry the fit is good (x 2 /dof= 873/880 = 0.99) and best 
fit parameters kT e = 190 ± 130 keV and r = 0.22^q'q5 are 
found. 

The measured iron Ka equivalent width for the 2012 
and 2013 observations (EW~35-200 eV, Tabic [2]) are un- 
usually large given the observed Compton reflection for the 
line to originate from reflection from Com pton thi ck ma- 
terial. We therefore confirm the result in Malaguti et al. 
( 1999) from the low flux BeppoSAX observa tion and with 
Suzaku ( Reeves et al ]| 20061 : iRlvers et al.ll2014l l. If the Comp- 
ton reflection and the iron Ka emission line arise from the 
same distant, Compton-thick material, they can be self- 
consisten tly modelled using the pexmon model in Xspec 
dNandra et al.lf2007l ). with the inclination angle fixed at 60 
degrees. This leads to a best fit x 2 /dof=930/884=1.06 with 
visible residuals around iron Ka. Leaving the iron abun- 
dance free to vary a significant improvement in the fit is 
found (Ax 2 = 34 for one additional variable parameter) 


and we find a Ap e > 22 with respect to the solar value. No 
significant variation between photon index, reflection frac- 
tion, absorbing column densities, cross-calibration factors 
and best fit parameters discussed above is found. Such a 
large iron abundance is unrealistic and therefore we discard 
this scenario. 

The alternative scenario we consider is the one in which 
the line is produced by Compton-thin material (Ah = 
10 22 — 10 23 cm” 2 ) that do es not contribute significantly 
to the Compton reflection (iMatt et al] 120031 ') . We model 
the Compton reflection and iron fluorescent l ines with MY- 
Torus scattered an d line components dMurphv fe Yaaoobl 
120091 ; lYaaoobl l2012h . We add a further component to re- 
produce absorbing material along the line of sight. The as- 
sumed geometry is a torus of gas and dust with a 60 degree 
opening angle. When we apply this model to the 2012 and 
2013 Swift+NuSTAR simultaneous data, we leave the nor- 
malization of the primary continua and the column den- 
sity along the line of sight as variable parameters. Nor- 
malizations and column densities of the scattered and line 
components are linked and free to vary. The fit is good 
(x 2 /dof= 907/884 = 1.02) and we find no variations from 
best fit parameters discussed above are found. The best fit 
value for the column density of the scattering material pro- 
ducing the iron Ka is Ah = 2.0 ± 1.1 x 10 23 cm” 2 . The 
ratios between the scattered and primary components’ nor- 
malizations is consistent with unity, the stand ard MyToru s 
configuration (coupled reprocessor model in lYaaooblboid h 
in which the torus is not aligned along the line of sight. 

The value of the column density of the iron K a emitti ng gas 
is con sistent with the estimate presented in IBianchi et alJ 
(2003) (A h — 3 x 10 23 cm” 2 ) for the case of NGC 7213, 
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Figure 4. 3-10 keV phenomenological best fit. EPIC-Pn data are 
in black, Suzaku data from the 2005 and 2013 observations are 
in red and green, respectively. NuSTAR FPMA and FPMB data 
from the 2012 observation are in light and dark blue while NuS- 
TAR FPMA and FPMB spectra obtained in 2013 are in orange 
and yellow, respectively. 


where the authors assumed a Broad Line Region covering 
factor f c = 0.35, an EW~ 100 eV and a photon index 
r = 1.69. This suggests that the two sources do indeed 
present similar features in the X-ray band. 
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3.3 Time history of the iron Ka line 

We apply the model described in Sect. 13.21 (an absorbed 
primary continuum and two Gaussian lines) to the archival 
3-10 keV data sets. This phenomenological fit is intended 
for studying the Fe Ka equivalent width and flux on time 
scales of months and years. 

We first fit data with high SNR from XMM-Newton, 
Suzaku and NuSTAR with the above model: the fit is good 
(x 2 /dof=628/544=1.15) and no strong residuals are present 
across the energy band (Fig. [4j. The best fit parameters are 
reported in Table 2. If we introduce a reflection component 
in our fit, no variation in the parameters of the lines is found 
with respect to the best fit values. 

Then, we fit the 3-10 keV spectra of the nine Swift snap- 
shots of the source, fixing the photon index to F = 1.7 (as 
inferred from the high energy data analysis). Further param- 
eters such as centroids and widths of the Fe Ka emission 
line are fixed to 6.4 keV and 60 eV (as inferred from the 
XMM and Suzaku analyses), respectively, where SNR. was 
too poor to accurately measure these parameters. Best fit 
X 2 /dof values are reported in Table 2. 

Fig. [5] shows the time evolution of the equivalent width 
of the Fe Ka line flux and the observed 3-10 keV flux of 
the source. The object has a variable intrinsic emission in 
the 3-10 keV energy range: we measure a factor of ~7.5 
between the 2003 XMM-Newton and 2012 NuSTAR obser- 
vations, while we find a factor ~ 1.5 between the Swift ob- 
servation on 2013-02-03 and the NuSTAR observation on 
2013-02-14 on a 10 day time scale. 

The variability of the Fe Ka emission line carries in- 
formation about the distance and ionization state of the 
emitting material (width and energy centroid). Indeed, the 


Figure 5. Time evolution of iron Ka parameters and 3-10 keV 
observed flux of the source (not corrected for absorption) . Equiv- 
alent widths are in eV units, line fluxes are in 10 -5 ph cm~ 2 s -1 
units and observed 3-10 keV fluxes in 10” 11 erg cm 2 s 1 units. 
Solid and dashed horizontal lines represent mean and standard 
deviations, respectively. 


measured Fe Ka flux in the 2003 XMM observation is sig- 
nificantly lower than the one measured with Suzaku three 
years later. 

In Fig. [6] we plot Fe Ka fluxes and equivalent widths 
(EW) vs the observed 3-10 keV flux (F^Lm), respectively. 
We examined the correlations between EW and F|_ 10 and 
between the Fe Ka flux and F-f-io by performing a linear fit 
using only high SNR observations (i.e. removing Swift data 
points) . 

We find a best fit relation EW = (—8.2 ± 1.6) x 
F 3 _ 10 + (163 ± 20) with a Spearman’s rank correlation co- 
efficient p = —0.88 ± 0.06 and a null-hypothesis-probability 
of 5.5 x 10~ 3 . The other best fit relation is Fx a = (0.42 ± 
0.17) x F 3 _ 10 + (5.2±1.0) with a Spearman’s rank correlation 
coefficient p = 0.67 ± 0.09 and a null-hypothesis-probability 
of 9.1 x 10“ 2 . Both best fit curves are shown in Fig. [6] When 
we introduce the Swift data points in the fit, we obtain 
less statistically significant relations, with null-hypothesis- 
probabilities of 3.9 x 10 -2 and 0.37 for the EW~F 3 _ 10 and 
Fko,-F§_ 10 relations, respectively. 
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Figure 6. Left panel: Iron K a flux (in 10“ 5 ph cm -2 s — 1 units) plotted against intrinsic 3-10 keV flux (in 10 11 erg cm -2 s — 1 units). 
Right panel: Iron Ka equivalent widths (in eV units) plotted against intrinsic 3-10 keV flux (in 10' 11 erg cm -2 s — 1 units). Dashed black 
lines represent best fit relations, when a linear fit is applied using only high SNR observations (i.e. removing Swift data points). 


4 DISCUSSION 


4.1 Properties of the hot corona 

Recently, the NuSTAR’s high sensitivity above 10 keV 
has led to high-energy cutoff measurements in a num- 
ber of nearby Seyfert galaxies: IC 4329A (178+40 keV, 
iBrenneman et all 120141). S WIFT J2127.4+5654 (108tiJ 


keV, 


Marinucci et, al.l feoidl h Ark 120 ( E c > 190 keV, 


iMat't et alj 20141) and 3C 382 (214+gg 7 keV, Ballantyne et 


al., submitted). The broad band analysis of NGC 2110 pre- 
sented in Sect. nm leads to a lower limit on the high-energy 
cutoff, E c > 210 keV. The unabsorbed 2-10 keV luminosity 
ranges between Lx =0.4-3.5x 10 43 erg s -1 (considering the 
2003 XMM and 2012 NuSTAR observations as low and high 
flux states, respectively ). If we use the 2-10 keV bolometric 
corrections presented in lMarconi et al] d2004l ). the bolomet- 
ric luminosity ranges between Lbol =0.6-9xl0 44 erg s ~ 4 


A black hole mass of Mbh — 2 x 10 8 Mq ll Merloni et al.l 
120031 ; iMoran et al.|l2007i ) leads to an Eddington luminosity 
LEdd — 2.4 x 10 46 erg s _1 . Therefore, for NGC 2110 we es- 
timate I,boi/TEdd — 0.25-3.7 x 10 -2 . This value is at the 
lower end of the distribution recently measured in CAIXA 
(Catalog of AGN in the XMM-Newton Archive) for a sample 
of 156 AGN dBianchi et alJl2009l) . 

NGC 2110 is the second Seyfert galaxy kwown to unam- 
biguously lack reprocessed emission from distant Compton- 
thick material and to show an iron Ka emission line likely 
produced by Compton-thin material. NGC 7213 is the 
other sourc e known to have simil ar properties l|Bianchi et aid 
120031 . l2008l : iLobban et al.1 120101 ) . The latter is also accret- 
ing at a low Eddington rate (I/boi/L/Edd ~ 3 x 10 -3 ) and 
IBianchi et al.1 d2008l ) reported a Broad Line Region origin 
for the iron Ka line. Interestingly, even though the two 
sources show remarkably similar spectral features in the X- 


rays, a high energy cutoff E c = 95 + 2 Q keV was detected 
in NGC 7213 in the simultaneous XMM-BeppoSAX data 
dBianchi et al.ll2004l ). 

In the near future, more AGN will be observed by NuS- 
TAR and we will be able to investigate the coronal proper- 
ties (temperature, geometry, link to the reflected emission 
from the accretion disk) with greater detail in additional low 
accretion rate objects. 


4.2 Iron Ka temporal behavior 

We analyzed the 3-10 keV spectra of NGC 2110 from sev- 
eral X-ray observatories, spanning a period of 16 years. 
Large variations in the observed flux of the source and in 
the Fe Ka line flux and equivalent width are apparent. If 
the Fe Ka line is produced by cold, distant matter we ex- 
pect a constant line flux and an EW linearly anti-correlated 
with the intrinsic flux of the source. On the other hand, 
if the Fe Ka line-emitting material is closer than distances 
corresponding to the time scales between the observations 
we expect a constant EW and a line flux linearly corre- 
lated with the illuminating continuum. We instead find a 
correlation between the Fe Ka flux and the 3-10 keV ob- 
served flux Fxa oc 0.42 x (E 3 L 10 ) and an anti-correlation 
between the Fe Ka EW and the 3-10 keV observed flux 
EW oc —8.2 x (F 3 _ 10 ). We therefore propose a scenario 
where the Fe Ka line is the sum of two distinct components, 
one constant and produced from material distant from the 
nucleus (the putative “torus”) and the other one variable 
and linearly correlated with the primary flux that is likely 
associated with the Broad Line Region. In this scenario the 
intercept of the relation between the Fe Ka flux and the 3-10 
keV observed flux gives the amplitude of the constant com- 
ponent. We find a value of (5.2 ± 1.0) x 10 ~ 5 ph cm -2 s _1 , 
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perfectly consistent with the resolved Chandra component 
(Table [2J . The Compton reflection fraction associated to 
this constant iron Ka component is consistent with the one 
found in the NuSTAR T Swift analysis fSect. lT2l) . 

The emitting structure responsible for the variable com- 
ponent of the Fe Ka line, clos er to the nu c leus t han the 
constant one, was discussed in I Evans et al.l (120071 '). These 
authors detected a modest broadening of the Fe Ka line 
in multi-epoch coadded HEG spectra, with a line width 
a — 4500^2200 km s _1 . We speculate that this struc- 
ture could be the same responsible for the broad, double- 
peaked Ha lines (FWHM ~ 13,000-17,000 km s -1 ) de- 
tected in optical spectro-polarimetri c observations of this 
source jMoran et al.ll2007l ; I li arJlioTol ). 

The Astro-H satellite, with its unprecedented combina- 
tion of spectral resolution and collecting area at 6-7 keV, 
will allow us to resolve the iron line profile at 6.4 keV in 
NGC 2110 and determine the location of the circumnuclear 
emitting material and its Fe abundance. 


4.3 A closer look to the circumnuclear 
environment 

It is interesting to note that NGC 2110 appears to show 
features at infrared wavelenghts arising from complex cir- 
cumnuclear regions. If we consider the mid-in frared lumi- 
nosity L\ 2 fim = 1.0 X 10 43 erg/s, reported in IHonig et al.l 
(|2010f ) and the 2-10 keV luminosity range discussed above, 
NGC 211 0 lies on th e Lmir oc (Lx ) 1 ' 11 relation inferred in 
iGandhi et al.l d2009l h who analized a sample of 42 Seyfert 
galaxies. This indicates that the geometry of the infrared 
emitter in NGC 2110 does not differ dramatically from other 
Seyfert galaxies. To produce enough IR continuum emission, 
the main requirement is enough dust with a sufficient cover- 
ing factor to reprocess the intrinsic AGN power. This could 
be satisfied by a number of models, incl uding a Con r pton- 
thin torus, or an extended dusty wind (IHonig et al.l l2()fid . 
and references therein). In this way, one could potentially 
have strong dust emission, with very little accompanying 
X-ray Compton reflection. 

I 11 addition, NGC 2110 is one of the few Seyfert galaxies 
(together with N GC 721 3 ) that also shows Silicate dust fea- 
tures in emission dHonig et al.l I2OI0I ') . The origin of Silicate 
emission is not fully understood even in Type 1 AGN, with 
non-standard dust grain properties, emission from a clumpy 
torus with a relatively small number of dust clumps, or 
emission from an extended dusty component i n the Narrow 
Line Region, all invoked as possib l e sour ces dSturm et aid 
120051: IShi et aid [200(3) . iMason et al.l d2009lf discussed in de- 
tail such a feature and found that the mid-IR component 
in NGC 2110 cannot be extended more than 32 pc, rul- 
ing out extended refl ecting clouds as seen in NGC 4945 
dMarinucci et al.ll20L3h We conclude that a standard dusty 
torus/dusty wind model, but with a gas column density in 
the Conrpton-thin regime could explain the X-ray and mid- 
IR characteristics of NGC 2110. 

iBeckmann fe Do Caol (|2010h reconstructed the spectral 
energy distribution of NGC 2110 using simultaneous IN- 
TEGRAL and Swift data taken in 2008 and 2009, report- 
ing features usually shown by radio-loud sources. However, 
we do not confirm the flat photon index and cutoff ener- 
gies they report. We use the radio fluxes at 6 cm (Fecm — 


165 mjy, Gnffithfrtal 


119951) and 20 cm (F 20 cm — 300 


nrjy, Brown et al. 2011 ) and compare them with the 2- 
10 keV fluxes we found. No significant deviations from re- 
lations between radio and X-ray em ission usually found 
in radio-quiet Seyf erts are present (IPanessa et_ah| _ 2007 : 
iBiemchi et al.l I2009L and references therein - ). lEvans et al.l 
( 200611 analyzed the Chandra , HST, and VLA imaging ob- 
servations and found a small radio jet (extended by ~ 5" 
accross the nucleus). However, the authors discarded the 
possibility of a synchrotron origin for the X-ray emission 
in NGC 2110, since the radio and X-ray emission are not 
spatially coincident. 


5 SUMMARY AND CONCLUSIONS 

We report a multi-epoch X-ray spectral analysis of the bright 
Seyfert 2 galaxy NGC 2110, spanning a period of 16 years. 
We focus on recent observations of the source with NuSTAR 
in 2012, when the source was at the highest flux level ever 
observed, and in 2013, when the source had more typical 
flux levels. 

Our results can be summarized as follows: 

• a high energy cutoff E c > 210 keV has been inferred, 
with an upper limit on the Compton reflection contribu- 
tion of R < 0.14, confirming results f rom past high-energy 
BeppoSAX and Suzaku observations IlMalaguti et al.l 1 19991 : 
iReeves et al.ll2006i : [Rivers et al. I mi); 

• when multi-epoch data are considered, we find a cor- 
relation between Fe Ka EW and intrinsic 3-10 keV flux 
(EW oc —8.2 x F 3 - 10 ) and an anti-correlation between Fe 
Ka flux and intrinsic 3-10 keV flux ( Fk& oc 0.42 x F3L10); 

• the Fe Ka line is likely the sum of two components: 
one constant (originating from distant Conrpton-thick ma- 
terial) and the other one variable and linearly correlated 
with the source flux (from matter at distances compati- 
ble with the Broad Line Region). Using MyTorus self- 
consistent modeling, we find that the line could be produced 
by scattering material with a global covering factor of 0.5 
with Ah = 2.0 ± 1.1 x 10 23 cm -2 ; 

• the source presents remarkably similar features to the 
low accretion rate Seyfert 1 galaxy NGC 7213 in the X- 
ray band (lack of Compton reflection, contribution from 
Compton-thin material to the Fe Ka line emission) and in 
the infrared , where Silicate dust emission was reported in 
iHonig et al. | ( 20 1 ni l. 
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